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ABSTRAGT 

The [G II] 158 /im fine-structure line is known to trace regions of active star formation and is the main 
coolant of the cold, neutral atomic medium. In this Letter, we report a strong detection of the [C II] 
line in the host galaxy of the brightest quasar known at z > 6.5, the Pan-STARRSl selected quasar 
PSO J036.5078-1-03.0498 (hereafter P036-I-03), using the IRAM NOEMA millimeter interferometer. Its 
[C II] and total far-infrared luminosities are (5.8±0.7) x 10® Lq and (7.6±1.5) x 10^® Lq, respectively. 
This results in a L[cii]/Atir ratio of ^ 0.8 x 10“®, which is at the high end for those found for active 
galaxies, though it is lower than the average found in typical main sequence galaxies at z ^ 0. We 
also report a tentative additional line which we identify as a blended emission from the 822 — 813 and 
523 ~ 432 H2O transitions. If confirmed, this would be the most distant detection of water emission to 
date. P036-I-03 rivals the current prototypical luminous J1148-I-5251 quasar at z = 6.42, in both rest- 
frame UV and [C II] luminosities. Given its brightness and because it is visible from both hemispheres 
(unlike J1148-I-5251), P036-I-03 has the potential of becoming an important laboratory for the study 
of star formation and of the interstellar medium only ~ 800 Myr after the Big Bang. 

Subject headings: cosmology: observations — quasars: emission lines — quasars: general 


1. INTRODUCTION 

The study of star formation in the first galaxies back 
in the epoch of reionization (z > 6.5) is one of the main 
challenges in current observational cosmology. 

Several groups have tried to identify the [C II] 158 A^m 
fine-structure transition line (hereafter [C II]) in the first 
spectroscopically identified galaxies at these redshifts. 
The [G II] line is the dominant coolant of interstellar 
neutral gas and is one of the brightest lines in the spec¬ 
trum of star-forming galaxies, accounting for up to ~ 1% 
of the total far-infrared (FIR) luminosity in z ~ 0 main 
sequence galaxies. Studies of the [C II] line may also 
provide key insights into galaxy kinematics at the high¬ 
est redshifts (see Garilli & Walter 2013, for a review). 
However, thus far, there are no convincing detections 
of [C II] emission from star-forming galaxies during the 
epoch of reionization (Walter et al. 2012; Ouchi et al. 
2013; Kanekar et al. 2013; Gonzalez-Lopez et al. 2014; 
Ota et al. 2014; Schaerer et al. 2015; Watson et al. 2015; 
but see also the recent work of Maiolino et al. 2015). 

On the other hand, the [G II] line has been identified in 
a number of more extreme objects such as submillimeter 
galaxies and quasar host galaxies at z ^ 6 (Maiolino 
et al. 2005; Walter et al. 2009a; Riechers 2013; Wang 
et al. 2013; Willott et al. 2013, 2015). These objects are 
most likely the progenitors of massive early-type galaxies 
seen in the present universe. Thus, luminous quasars are 
thought to be important tools to pinpoint the locations 
of these extreme galaxies in the early universe. Indeed, 
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the [G II] line has been detected even in the host galaxy 
of tbe currently most distant quasar known at z = 7.1 
(Venemans et al. 2012). 

For more than a decade, the quasar SDSS J1148-I-5251 
at z = 6.42 (hereafter J1148-I-5251), discovered by Fan 
et al. (2003), has been the brightest quasar known at 
z > 6 (Mi 45 o = —27.8). As such, it is by far the best 
studied quasar at these redshifts, being fundamental in 
shaping our current understanding of the universe when 
it was less than 1 Gyr old. It was the first source de¬ 
tected in [G II] at z > 0.1 (Maiolino et al. 2005), and for 
several years it remained the only z > 6 quasar detected 
also in GO emission (Walter et al. 2003, 2004; Bertoldi 
et al. 2003b; Riechers et al. 2009). The [G II] line in 
this quasar has been studied in detail (e.g., Walter et al. 
2009a; Maiolino et al. 2012). Neutral carbon [C I] has 
also been detected (Riechers et al. 2009), while deep lim¬ 
its on the [N II] line could also be placed (Walter et al. 
2009b). 

Recently, and after more than 10 years, three quasars 
that rival (and in one case even surpass) J1148-1-5251 in 
UV luminosity have been discovered. These quasars are 
SDSS J0100-k2802 (z = 6.30; M^so = -29.3; Wu et al. 
2015), ATLAS J025.6821-33.4627 (z = 6.31;Mi45o = 
-27.8; Carnall et al. 2015), and PSO J036.5078-k03.0498 
(hereafter P036-I-03; z = 6.527; M1450 = —27.4; Vene¬ 
mans et al. 2015). 

In this Letter, we focus on P036-I-03, which is the 
brightest quasar known at z > 6.5 (see Figure 3 in 
Venemans et al. 2015). We report a bright detection 
of the [G II] line with the IRAM NOEMA interferome¬ 
ter. We also detect the underlying dust far-infrared con¬ 
tinuum emission and present a tentative detection of a 
water transition in this quasar. As a result, P036-I-03 
(that is also accessible from telescopes in the southern 
hemisphere) competes with the current archetypal high- 
redshift quasar J1148-I-5251 in both rest-frame UV and 
[G II] luminosities. 
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We employ a cosmology with Hq = 67.8 km 
Mpc“^, = 0.307, and f^A = 0.691 (Planck Collab¬ 
oration et al. 2014). This implies that for z = 6.5412, 
the proper spatial scale is 5.6 kpcarcsec”^ and the age 
of the universe is 833 Myr, i.e., 6% of its present age. 

2. OBSERVATIONS AND RESULTS 

We use the IRAM NOEMA interferometer to observe 
the [C II] line in the quasar P036-f03, as a director’s 
discretionary time proposal (ID: E14AG). The tuning 
frequency is 252.496 GHz which corresponds to the ex¬ 
pected [C II] frequency {I'rest = 1900.539 GHz) based 
on the quasar redshift estimated from the Mgii line: 
-^Mgii = 6.527 (Venemans et al. 2015). The observations 
were carried out in 2015 February 09 and 10 in the D 
configuration with a total time on source of only 3.55 
hr (5 antenna equivalent). The synthesized beam size is 
2.2" X 1.6". The data were calibrated and analyzed with 
the IRAM GILDAS^ software package. 

The hnal NOEMA spectrum of P036-I-03 is shown in 
Figure 1. Two lines are visible, a bright one, and a 
second fainter line on top of the underlying dust con¬ 
tinuum. We ht two Gaussians plus continuum to the 
data and identify [C II] as the brightest of these lines. 
The [G II] line peaks at 252.02 ± 0.06 GHz, yielding 
a redshift of Z[cii] = 6.54122 ± 0.0018 (i.e., redshifted 
by 566 ± 72kms“^ with respect to the Mgii redshift). 
The line has a FWHM= 360 ± 50kms“^ and a veloc¬ 
ity integrated flux of 5.2 ± 0.6Jykms“^. This corre¬ 
sponds to a [G II] luminosity of T[cii] = (5.8 ± 0.7) x 
10® Lq which is a factor of ~1.4 brighter than the [G II] 
line in J1148-I-5251 (4.2 ± 0.4 x 10® L©; Maiolino et al. 
2005; Walter et al. 2009a). The second line peaks at 
254.45±0.14 GHz, i.e., blueshifted by —2863±178kms“^ 
with respect to the [C II] line. This fainter line has a 
FWHM= 160 ± 100 kms“^ and an integrated flux of 
1.2 ± 0.7Jykms“^. At our resolution, the most likely 
identification for this line is a blended water vapor emis¬ 
sion from the adjacent 822 — 813 and 623 — 432 H 2 O tran¬ 
sitions {vrest = 1919.359 GHz and Vrest = 1918.485 GHz, 
respectively). This tentative blended water emission has 
a velocity offset of 576 ± 165 kms“^ with respect to its 
predicted value by the Mg II redshift. This offset is in 
good agreement with the velocity offset found for the 
[G II] line, reinforcing the identification of the water 
emission. Figure 2 shows the maps of the far-infrared 
continuum emission and the continuum subtracted [G II] 
emission. 

Fitting the continuum flux density yields a value of 
2.5 ± 0.5 mjy. We convert the measured dust contin¬ 
uum flux density to far-infrared luminosity, Tfir, by 
modeling the FIR emission with an optically thin gray 
body spectrum and scale the model to match the ob¬ 
served dust continuum flux density. We use typical pa¬ 
rameters in the model to reproduce the SED for high- 
redshift quasar host galaxies, i.e., a dust temperature 
of Td = 47 K and a emissivity index of /3 = 1.6 (fol¬ 
lowing for example, Venemans et al. 2012; Willott et al. 
2015). Integrating the model from rest-frame 42.5/rm 
to 122.5/rm we obtain the following LpiR for P036-1-03: 
Afir = (5.4 ± 1.1) X 10^® L©. Integrating the model 
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from rest-frame 8 /rm to 1000 fim we get a total infrared 
luminosity, Atir, of Atir = (7.6 ± 1.5) x 10 ^® A©. We 
can use the total infrared luminosity to estimate the star 
formation rate using the relation SFRtir(M 0 yr“^) = 
6mf X 1.0 X 10-1® Atir (A©) (Garilli & Walter 2013, 
and references therein), where Smf depends on the stel¬ 
lar population. We use 6mf = 1.0 which is appropri¬ 
ate for a Chabrier IMF. For P036-I-03, this relation re¬ 
sults in SFRtir = 760 ± I 5 OM 0 yr“^. Alternatively, 
Sargsyan et al. (2014) show that even in sources domi¬ 
nated by AGN, the star formation rate can be estimated 
from A[cii] to a precision of ~ 50% with the relation 
SFR[cii](M 0 yr“^) = IO'^Ajqh] (A©). This relation 
yields SFRjcn] = 580 ± 7 OM 0 yr“^, broadly consistent 
with the estimate above. 
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Fig. 1.— NOEMA spectrum at 1mm (observed) wavelengths 
of P036+03. The tuning frequency was centered on the expected 
position of [C II] based on the redshift determined from the Mgii 
line (uobs = 252.496 GHz). This tuning frequency is shown here as 
zero point of the velocity scale. The red solid curve is a fit to the 
data consisting of two Gaussians plus continuum. The [G II] line is 
detected at high significance and is redshifted by 566 di 72kms~^ 
with respect to the Mgii redshift = 6.54122, ^Mgii = 6.527). 

A second, fainter line is observed at 254.45 di 0.14 GHz. We find 
that the most likely identification for this line is a blended emission 
from the adjacent 322 — 3i3 and 523 ~ ^32 H2O emission lines (see 
text). 

3. DISCUSSION 

3.1. [CII] - IR luminosity relation 

It has been shown that in local star-forming galax¬ 
ies the ratio between the [C II] and IR luminosities 
is about 0 . 1 % - 1 % but it decreases in galaxies with 
larger IR luminosities (Ltir > 10^^ Lq) such as LIRGs 
and ULIRGs. This has been known as the “[G II] 
deficit”. Similar deficits have been observed in other 
far-infrared lines such as [O I] 63.2 ^m, [O I] 145 /.im, and 
[N II] 122 yiva. As a consequence, these trends are now be¬ 
ing referred to as the “far-infrared lines deficit” (Gracia- 
Garpio et al. 2011). In Figure 3 we show the ratio of 
















































[C II] and dust in a bright quasar at z = 6.54 
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Fig. 2.— NOEMA cleaned maps of the 2: = 6.5412 quasar P036+03. The red cross in each panel shows the near-infrared location of 
the quasar. The beam size of 2.2" X 1.6” is shown at the bottom left of each panel. Left: Far-infrared continuum emission obtained 
from line-free channels in the frequency range: 254.274GHz - 252.524GHz (i.e., 2100kms“^, see Figure 1). The contours correspond to 
—3cr, —2cr, —1(7 (white dashed lines), Icr, 2t7, 3(7, 4a, and 5a (black solid lines), with a being the rms noise of 0.418mjy beam“^. Right: 
Continuum subtracted [C II] emission integrated from the channels 251.372GHz to 251.782 GHz (i.e., 750kms“^). The contours are —3 (t, 
—2a, —la (white dashed lines), la, 3a, 5a, 7a, and 9a (black solid lines), with a = 0.85mjy beam“^. The excess of 3a sources in the 
maps is due to the modest u, v coverage in our data. 


[C II] to TIR luminosity for different type of sources. 
The “[C II] deficit” can be clearly appreciated in the lo¬ 
cal samples: star-forming galaxies (Malhotra et al. 2001), 
LIRGs (Diaz-Santos et al. 2013), and ULIRGs (Luhman 
et al. 2003). The sample of IR-luminous star-forming 
galaxies, submillimeter galaxies, and quasars at z > 1 
compiled by De Looze et al. (2014) and Brisbin et al. 
(2015) shows a larger scatter in the [G II] to TIR lumi¬ 
nosity ratio (^ 10“^ — 10“^'®). The sample of z ~ 6 
quasars with high IR luminosities studied by Maiolino 
et al. (2005); Walter et al. (2009b); Wang et al. (2013) 
shows small [C II] to FIR luminosity ratios, similar to 
what is observed in local ULIRGs (~ 10“^'®). One 
interpretation of such low ratios could be in part due 
to the central AGN contributing to the IR luminosity 
(Wang et al. 2013), but note that there is a big liter¬ 
ature dealing with possible explanations for the deficit 
(e.g., Malhotra et al. 2001; Gracia-Carpio et al. 2011; 
Dfaz-Santos et al. 2013). Willott et al. (2015) studied a 
sample of z ^ 6 quasars two order of magnitudes fainter 
in the infrared (Ttir = 10^^“^^ Lq). These quasars 
show ratios consistent with the local star-forming galax¬ 
ies (~ 10“^ —10“^'®). Venemans et al. (2012) report that 
the quasar J1120-I-0641 at z = 7.1 has a TIR luminos¬ 
ity that is between the Willott et al. (2015) sample and 
the IR-bright z ^ 6 quasars. The L[Cii]/Atir ratio for 
J1120-1-0641 is more consistent with local star-forming 
galaxies than to ULIRGs. The quasar of the present 
study, P036-I-03, has a Ltir comparable to other z ~ 6 
quasars (Ltir 10^^ L©). However, it has a signih- 
cant larger L[ciii/.bTiR ratio (^ 0.8 x 10“^), similar to 
J1120-I-064I, and consistent with the lower-end of local 
star-forming galaxies. This ratio is 2.7 times larger than 
the ratio reported for JI148-I-5251. 

Together, this emphasizes the differences between the 
host galaxies of the most distant quasars. However, the 
number of quasars with [C II] and dust far-infrared con¬ 
tinuum information is still to small to derive statistically 


meaningful trends or correlations. 
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Fig. 3.— The ratio of [C II] to total IR luminosity (TIR; 8/^m 
to 1000/im) vs TIR luminosity. The sources at 1.0 < z < 6.4 in¬ 
clude the compilation of star-forming galaxies, submillimiter galax¬ 
ies, and quasars by De Looze et al. (2014) and 8 star-forming 
galaxies recently reported by Brisbin et al. (2015). The 5.7 < 
2 < 6.5 quasars are from Bertoldi et al. (2003a); Maiolino et al. 
(2005); Walter et al. (2009a); Wang et al. (2013); Willott et al. 
(2013, 2015). The z > 6.5 quasars are P036+03 (this work) and 
J1120-h0641 at 2 = 7.1 (Venemans et al. 2012). Luminosities have 
been calculated with the cosmology used in this paper. The TIR 
luminosities for the 2 > 5.7 quasars have been consistently cal¬ 
culated as for P036-I-03 (see text Section 2). The FIR (42.5/^m 
to 122.5 fim) luminosities from Diaz-Santos et al. (2013) have been 
converted to TIR luminosities using a common conversion factor of 
1.75 (Calzetti et al. 2000). Error bars are only plotted for P036-I-03 
to enhance the clarity of the figure. 


3.2. Water vapor detection 

The water vapor lines arise in the warm, dense phase of 
the ISM. They are one of the main coolants of this phase, 
thus playing a key role in the fragmentation and collapse 
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of gas clouds (Carilli & Walter 2013). To date, there are 
only a handful of water vapor lines reported at high red- 
shift, mostly in strongly magnified sources (e.g., van der 
Werf et al. 2011; Omont et al. 2013) and in one extreme 
submillimeter galaxy at z = 6.3 (Riechers 2013), being 
the most distant detection of interstellar water to date. 
Recently, Bialy et al. (2015) show that high abundances 
of water vapor could exist in galaxies with extremely low 
metallicity 10 “^ solar), as expected for the first gen¬ 
erations of galaxies at high redshifts. 

As discussed in Section 2, we identify a tentative 
detection of blended water vapor emission from the 
322 — 3 i 3 and 823 — 432 H 2 O transitions (see Figure 
1). Owing to the low significance of the detection is 
difficult to make any strong conclusions, especially be¬ 
cause the water emission is complex and hard to in¬ 
terpret if only one transition is measured (Gonzalez- 
Alfonso et al. 2014). However, additional water tran¬ 
sitions could be detectable in this source. For example, 
the 321 — 3 i 2 H 2 O transition (i^rest = 1162.912 GHz) at 
2 ; = 6.5412 is shifted into the 2 mm band and could si¬ 
multaneously be observed with the GO(J= 10 —9) line 
(i'rest = 1151.985GHz). The detection of these transi¬ 
tions would lead to an improved excitation model for one 
of the earliest galaxies known in the universe and to the 
confirmation of the most distant water detection to date. 

4. SUMMARY 

We detect a bright [C II] emission and the underly¬ 
ing continuum emission in P036-I-03, the most luminous 
quasar known so far at 2 : > 6.5. A second line is ten¬ 
tatively detected which we identify with blended wa¬ 
ter emission from the adjacent 622 — 3 i 3 and 523 ~ 432 
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H 2 O transitions. If confirmed, this would be the highest- 
redshift detection of water vapor available to date. 

The [G II] to TIR luminosity ratio in P036-I-03 is ^ 
0.8 X 10“^, i.e., consistent with Iow-Ltih star forming 
galaxies and high compared to other IR-luminous quasars 
at z > 6 . This ratio is a factor ~ 3 higher than in 
J1148-I-5251, and consistent with the ratio in the most 
distant quasar known to date (Figure 3), emphasizing the 
diversity of quasar host galaxies in the early universe. 

P036-I-03 rivals the current archetypal high-redshift 
quasar, J1148-I-5251 (z = 6.42), in both UV and [C II] 
luminosities. Thanks to its brightness and given its con¬ 
venient equatorial location, P036-I-03 can be studied from 
facilities in both hemispheres (in contrast to the north¬ 
ern source J1148-I-5251). For example, P036-I-03 would 
be an ideal target for high-resolution ALMA imaging to 
study the ISM physics and to search for signs of mergers, 
outflows, and feedback in one of the earliest galaxies in 
the universe. This quasar will thus likely become one of 
the key targets to study star formation and the ISM in 
the early universe. 
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